Copyright © 1979 American Telephone and Telegraph Company 
The Bell System Technical Journal 
Vol. 58. No. 5, May-June 1979 
Printed in U.S.A. 



Digital Radio Outage Due to Selective 

Fading — Observation vs Prediction From 

Laboratory Simulation 

By C. W. LUNDGREN and W. D. RUMMLER 
(Manuscript received July 7, 1978) 

A statistical model (introduced in a companion paper) of fading on 
a radio path is used with laboratory measurements on a digital radio 
system to estimate the outage due to multipath fading, where outage 
is the time that the bit error rate (ber) exceeds a threshold. Over the 
range of ber of interest (10~ 6 to 1G~ 3 ), the calculated outage agrees 
favorably with the outage observed during the period for which the 
fading model was developed. It is further shown that the calculated 
outage, when scaled to a heavy fading month on the basis of single- 
frequency, time-faded statistics, agrees equally well with the outage 
observed on the same path during a heavy fading month. The agree- 
ment between measured and predicted outage substantiates the se- 
lective fading model. The prescribed laboratory measurements char- 
acterize the sensitivity of the radio system to selective fading. Thus, 
the methodology provides a useful basis for comparing the outage of 
alternative realizations of digital radio systems. 

I. INTRODUCTION 

Present interest in using high-speed common carrier digital radio 1-5 
has precipitated a need for estimating the performance of such systems 
during periods of selective (multipath) fading. This paper describes a 
method of characterizing a digital radio system in the laboratory which 
allows the outage to be accurately estimated. For a digital radio 
system, outage requirements are stated in terms of the number of 
seconds in a time period (usually a heavy fading month) during which 
the bit error rate (ber) may exceed a specified level; typically, 10" 3 or 
10~ 4 is appropriate to voice circuit application. 

The method is based upon a statistical channel model 6 developed 
from measurements on an unprotected 26.4-mile hop in the 6-GHz 
band in Palmetto, Georgia in 1977 using a general trade 8-psk digital 
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radio system as a channel measuring probe. The modeled fading 
occurrences were scaled to the basis of a heavy fading month using the 
occurrence of time faded below a level at a single frequency as the 
means of calibration. The bit error rate performance of the digital 
radio system was measured during the time period used for channel 
modeling and for an extended period corresponding to a heavy fading 
month. This same radio system was later subjected to a measurement 
program in the laboratory using a multipath simulator which provides 
a circuit realization of the fading model. The measured results are used 
with the channel model to determine the occurrence of channel con- 
ditions which will cause the ber to exceed a given threshold. Compar- 
isons on the basis of the modeling period and a heavy fading month 
show good agreement between calculated and observed outages for 
bers between 10" 6 and 10~ 3 . 

The properties of the fixed-delay channel model are reviewed briefly 
in Section II as a basis for describing the measurements and for the 
subsequent outage calculations. This three-parameter channel model 
is realized in the laboratory by an if fade simulator. The simulator 
and its use in obtaining the necessary laboratory data are described in 
Section III. The procedures to be followed in calculating outage times 
for a given ber are described in Section IV. Calculated and observed 
outage times are compared in Section V. Conclusions are provided in 
Section VI. 

II. MODEL DESCRIPTION— METHODOLOGY 

It has been demonstrated 6 that the complex voltage transfer function 
of a line-of-sight microwave radio channel is well modeled by the 
function 

H(u) = a [1 - 6e- >( "- Uo)T ] (1) 

with t fixed. A 6-GHz channel (30-MHz bandwidth) has been charac- 
terized statistically by the model with t = 6.3 ns. Such a channel has 
a power transfer function given by 

| H(u>) | 2 = a\l + b 2 -2b cos(oj - <o„)t] (2) 

and an envelope delay distortion function, i.e., the derivative of the 
phase of H (co) with respect to o>, given by 

6t(cos(w - o>o)t - b) 
Di03) = 1 + b* - 26 cos(u - a*)/ (3> 

In the following paragraphs, we summarize the properties of the model, 
the statistics of the model parameters, and the measurement objec- 
tives. 
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2.1 Fixed delay model 

A plot of the attenuation produced by the fixed delay model of eq. 
(1) is shown in Fig. 1. Since the delay is fixed at 6.3 ns, the spacing 
between nulls of the response, 158.4 MHz, is much larger than the 
channel bandwidth. The parameters a and b control the depth and 
shape of the simulated fade, respectively. The parameter f (=uo/27r) 
determines the position of the fade minimum or notch. Both the notch 
frequency, f Q , and the response frequency, f, are measured from the 
center of the 30-MHz channel for convenience. 

The model function of eq. (1) may be interpreted as the response of 
a channel which provides a direct transmission path with amplitude a 
and a second path providing a relative amplitude b at a delay of 6.3 ns 
and with a phase of coot + v (independently controllable) at the center 
frequency of the channel. This interpretation is represented in Fig. 2 
by a phasor diagram at to = 0, the center frequency of the channel. 
Varying the frequency, to, over the channel bandwidth (30 MHz) moves 
the angle of the interfering ray through an arc of about 60 degrees (2w 
X 30 MHz X 6.3 ns » tt/3), centered at the position shown. This 
diagram is useful for understanding the fade simulation; it also provides 
an alternate means of describing the position of the notch. The notch 
position may be specified by its frequency, f , or by <f>, the angle of the 
interfering path at the center of the channel. 




-100 



Fig. 1— Attenuation of channel model function, H(w) = a[l - b exp(-j(u - wo)t)], for 
t = 6.3 ns, a = 0.1, b = 0.7. 
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Fig. 2— Phasor diagram of a modeled fade, </> = 45°, b = 0.7. 

From Figs. 1 and 2, or eq. (1), it may be seen that varying a changes 
the overall level and varying b changes the shapeliness of the modeled 
fade. Furthermore, if the minimum is within the 30-MHz channel 
bandwidth ( |<J>| < 30°), the fixed delay model can generate notches 
with a wide range of levels and notch widths. With the minimum out 
of band, it can generate a wide range of combinations of levels, slopes, 
and curvatures within the channel bandwidth. 

2.2 Model statistics 

The statistics of the model parameters were obtained from a selected 
data base during which heavy fading activity was observed. 6 The 
distribution of 6 is best described in terms of B = -20 log(l - b). 
Figure 3 shows the distribution of B and the least-squares straight line 
fit to the distribution over the region where it best represents selective 
fading— between B values of 3 and 23 dB. The channel is described by 
B greater than 23 dB for less than 0.15 percent of the observed time 
which makes the distribution less certain beyond this point. At the 
other extreme, during the periods of time when there is little or no 
selective fading, the channel is characterized by values of B less than 
3 dB. Thus, the fitted line represents a lower bound on the distribution 
for B less than 3 dB. Since the fitted line has an intercept of 5400 
seconds, we may model the fraction of 5400 seconds during which B 
exceeded a value X by the probability distribution 

P(B>X) = e- x/38 . (4) 

Thus the probability of finding a value of B between X and X + dX is 

PB (X)dX = ^e- x/ ™. (5) 

The distribution of a is lognormal with a standard deviation of 5 dB 
and a mean that is dependent on B (or b). Hence, the probability that 
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Fig. 3 — Distribution of B for model data base period. 



A = —20 log a has a value between Y and Y + dY is given by 

dY 



p A (Y)dY = 



,-[y-A„(B)] 2 /50 



5V2tt 



(6) 



The relationship between Aq, the mean of the distribution, and B is 
given in Fig. 4. 

The distribution of f a is found to be independent of A and B. It is 
usually simpler to work with </> rather than f Q . The two variables are 
simply related in that </> is defined on the interval (—it, tt) and a 2.5- 
degree change in <j> corresponds to a 1.1-MHz change in f a . For the 
fixed delay model, the variable $ has been found to have a probability 
density function that can be described as uniform at two levels, with 
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Fig. 4 — Mean and standard deviation of the distribution of A = -20 log a as a 
function of B. 

values less than tt/2 being five times more likely than values greater 
than tt/2. Thus, we have the probability density function per degree 

as: 

|4>|<90 c 



P*(<t>) = 



1 
216 



(7) 



1080 



90° < I <f> I < 180. 



The functions in eqs. (5) to (7) can be used to determine the 
probability of finding a, b, and f in some region of a-b-f„ space. This 
probability can be converted to number of seconds in the observation 
period by multiplying by 5400 seconds. To convert this probability to 
the number of seconds in a month requires scaling the data base. The 
scaling may be obtained from Fig. 5 which shows, for several frequen- 
cies in the band, the time during the model data base period that the 
channel was faded below a given level. Distributions are shown for 
average power in the band and for power at selected frequencies at the 
center and near the edges of the radio channel. (Frequencies indicated 
are at if where the center frequency is at 70 MHz.) For the path used, 
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Fig. 5 — Amplitude distributions for model data base period. 

one expects the received voltage (single frequency) relative to midday 
average to be less than L for 72.5 X 10 4 L 2 seconds in a month. 7 For the 
data base used, the fading is best described by 48 X 10 4 L 2 hence, the 
data represent % of a fading month. To obtain outage on a seconds- 
per-heavy-fading-month basis, the probabilities calculated with eqs. 
(5) to (7) must be multiplied by 5400 X 1.5 or 8100 seconds. 

2.3 Outage estimation 

The fixed delay model described above can be simulated with an 
equivalent circuit laboratory measurement to determine the equip- 
ment response to multipath fading. Conceptually, one determines 
critical values of A and B for which a specified error rate is achieved 
for each fade notch position. In practice, it is difficult to maintain a 
constant ber; it is more expedient to fix b and vary the carrier-to-noise 
ratio (a) while plotting the ber. From the resulting curves, one may 
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compute critical contours of A and B for each prescribed notch location 
and ber. Using eqs. (5) and (6), the probability that A and B he on the 
high error rate side of a given critical contour may be calculated. 

By repeating this calculation for a uniform set of notch positions 
and using (7) to determine the probability weighting given to each and 
summing, one may estimate the probability of all selective fades that 
produce a ber exceeding the prescribed one. Multiplying this proba- 
bility by 5400 gives the outage time expected over the data base period; 
multiplying by 8100 gives the expected outage time per heavy fading 
month. 

The following section describes the laboratory measurement; Sec- 
tion IV describes the reduction of the measured curves and parameters 
to outage times. 

III. LABORATORY MEASUREMENTS 

Figure 6 illustrates stressing of a digital radio system by means of an 
if fade simulator. The simulator, which is inserted after linear if 
preamplification but before any high-gain amplification, shapes both 
the desired signal and the effective received noise. It is necessary to 
operate the simulator at a sufficiently large input carrier-to-noise ratio 
that the concomitantly shaped noise at its output remains a negligible 
contributor to degraded system performance throughout the operating 
range of interest. 

Within its restricted frequency range of operation, the if simulator 
is adjusted to achieve those specific shapes implied by Fig. 2. Although 
the measurements could have been made using an rf fade simulator, 
the choice of an if simulator was based primarily upon considerations 
of signal and noise levels, and the repeatability of adjustments. The 
following section describes an if shape-stressing measurement in the 
minimum detail necessary to qualify the data collected. 

3. 1 Representative IF two-path fade stressing measurement 

The block diagram of Fig. 7 illustrates an arrangement employing 
an if fade simulator and an if flat noise source. A pseudo-random test 
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Fig. 6— if fade stressing. 
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pattern modulates the 6-GHz radio transmitter whose output is nom- 
inally 5 watts ((g), in Fig. 7). The output spectrum is usually shaped by 
a bandpass filter following the microwave power amplifier to comply 
with FCC regulations. 

To enable back-to-back operation of the transmitter and receiver of 
a single repeater which normally operate on different radio channels, 
a radio test translator was employed. The translator output power was 
approximately —30 dBm (adjustable, at (S)) to simulate the unfaded 
received signal level (rsl) observed typically in the field. 

Assuming a linear rf-if conversion gain of 20 dB, the signal power 
at the input (©) to the if fade simulator is -10 dBm. The simulator 
incorporates low-noise linear amplification. A reference insertion loss 
for the main unfaded ray is 10 dB, including the output power summer. 
Hence the maximum desired signal power at the input to the main if 
amplifier (©) is -20 dBm. 

Assuming a 30-MHz receiver noise bandwidth and a current-art 
receiving system noise figure of 5 dB, the total system noise power is 
approximately —95 dBm, referred to the receiver's input port. This 
results in a flat receiver noise contribution of —85 dBm at input © to 
the main if amplifier. Consequently, the maximum attainable carrier 
to simulator-shaped rf noise ratio is 10 log (Co/N r f) = —20 — (—85) = 
65 dB. The noise contributed by the fade simulator amplifiers must 
not exceed —100 dBm, to be negligible. 

Flat if noise much larger than the unwanted and shaped system 
noise is added artificially at © and is adjusted in magnitude by a 
calibrated attenuator <3) to superpose thermal noise degradations upon 
the simulated selective fading degradations of the desired signal. One 
would ideally measure the added if noise power in the final predetec- 
tion bandwidth of the digital radio system, or twice the Nyquist 
bandwidth. It is more convenient in the laboratory to reference carrier- 
to-noise ratios to the output of the main if amplifier by using the - 
precalibrated agc voltage (assuming that wideband agc detection is 
employed), to measure both the unshaped signal and noise powers. 
The carrier-to-noise ratio at the detector would be higher — by the 
ratio of the system noise bandwidths that would be measured at the 
respective points. 

The noise source output in Fig. 7 may be adjusted so that an 
attenuator setting of dB (3) results in a noise power equivalent to 
that of the unfaded signal power (the attenuator is then calibrated 
directly in uncorrected Co/iV,/, in dB). As the if fade simulator is 
readjusted to achieve different prescribed fade shapes, its mean inser- 
tion loss may also change. The change in insertion loss is determined 
by monitoring the change in signal power at ©; the same loss increment 
(dB) must be added to the noise attenuator (3) to reestablish the dB 
reference. 
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3.2 IF two-path fade simulator 

Figure 8 illustrates splitting the desired if signal into an arbitrarily 
phased, adjustable "main" component and a "delayed" component 
fixed in delay (t ns) but adjustable in magnitude. The main component 
is further resolved into orthogonal components (inset to Fig. 8) using 
wideband networks exhibiting flat gain and well-behaved delay. A 
particular sum vector is constructed by adjusting the orthogonal 
components to establish a simulated fade notch frequency; in practice, 
the phase sense of 0- and 90-degree components are independently 
reversible, as indicated by the switches in the figure, for complete 
flexibility in notch frequency selection. 

The 6.3-ns fixed delay added to the delayed path imparts a phase 
shift of 159 degrees at the 70-MHz if center frequency. This is shown 
built out to 225 degrees, relative to the 0-degree transmission path, 
using a 66-degree wideband network of the same type. The delayed 
vector is fixed in direction opposite the midrange position of the 
adjustable main vector, corresponding to a channel-centered fade (<f> 
= degrees). 

Since 1/t = 158.4 MHz, a change of 1 degree in direction of the main 
vector corresponds to a frequency displacement of the fade notch 
location of 0.44 MHz. For <j> = —45 degrees, the notch location is 
displaced 19.8 MHz below the channel center (f = —19.8 MHz). The 
magnitude of the delayed component is then adjusted to achieve the 
desired notch depth. 
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3.3 Digital radio performance stressed by in-band selectivity and thermal 
noise 

The radio equipment was measured at uniformly spaced notch 
frequencies separated by 4.4 MHz (A<?> = 10 degrees). To fully charac- 
terize a period of variation in <f>, or f D , one would need to make 36 sets 
of measurements. Ideally, half may be omitted because of symmetry. 
For given values of A and B, the same error rate should obtain for a 
notch at a given frequency displacement above or below the channel 
band center. Variations in B ought not to have a significant effect for 
| <}> | greater than 90 degrees. It was determined that detailed measure- 
ments were required for nine different values of f a to characterize the 
digital radio tested. 

Using a wideband rf fade simulator in the field, the digital radio 
performance for out-of-channel notch locations was relatively inde- 
pendent of whether minimum or norimirurnum-phase fade simulations 
were employed. The nonminimum phase fade is modeled by eq. (1) 
with the sign of the phase term reversed. This leaves the amplitude 
[eq. (2)] unchanged, but reverses the sign of the envelope delay 
distortion [eq. (3)]. We conclude that the minimum phase channel 
model is sufficiently general for use in simulating the channel and in 
estimating performance. 

The if fade simulator was adjusted for each notch frequency, and 
the depth of notch was varied by adjusting the magnitude of the 
delayed component. Then various amounts of if thermal noise were 
added. Figure 9 typifies the performance data collected, bers are 
plotted versus the uncorrected if carrier-to-noise ratio (C/2V,/), for a 
constant fade notch offset from midchannel (f D = -19.8 MHz, r = 6.3 
ns). Each curve corresponds to a different notch depth (B = —20 
log(l-fc) dB), and hence a different amplitude and delay shape in the 
radio channel. Each curve is also identified with an in-band selectivity, 
defined as the difference between the maximum and minimum atten- 
uation present in the (25.3-MHz) channel bandwidth. The lower-left 
"baseline" curve presents the unshaped signal, flat fading performance 
obtained by adding only if thermal noise. This curve was verified 
(without the added if noise) by attenuating the received rf input 
signal in the back-to-back configuration. 

Consideration was given to matching the order of measurements to 
characteristics of the particular digital radio tested. For example, 
considerable scattering of data at low error rates can result from 
synchronizations involving different reference carrier phases. The au- 
thors elected to perform several synchronizations while observing the 
ber for each phase, and then chose that phase giving the worst 
performance.* Synchronization was accomplished at the low error rate 

* Because the phase information in the measured system was Gray coded and digital 
access was on a per-rail basis, one rail had twice the ber of the other two. All 
measurements in the field and in the laboratory were made on this worst-case rail. 
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Fig. 9 — High-speed digital radio if dispersive fade simulations, t = 6.3 ns, fa = —19.8 
MHz. 



(bottom) of each curve, and this phase relationship was maintained 
for all data points obtained for each curve. 

From the baseline curve of Fig. 9, a ber = 1 X 10" 6 obtains for 10 
\og(C/N if ) = 21.5 dB. For the digital radio system installed on the 
instrumented hop and reported in the figure, the measured flat fade 
margin for a threshold ber = 1 x 10 -6 was 40.5 dB This leads to an 
unfaded if carrier-to-noise ratio 10 \og(C /N if ) = 21.5 + 40.5 = 62 dB. 

From the baseline curve for a ber = 1 X 10" 7 , note that insertion of 
a fade whose notch depth is 6.5 dB results in four orders of magnitude 
degradation in ber performance; equivalently, an in-band selectivity 
of only 5.7 dB in 25.3 MHz results in a ber > 1 X 10" 3 . 

The asymptotic regions in Fig. 9, corresponding to high values of C/ 
Nif, are not normally presented in characterizations of this type; 
however, system outage depends primarily upon the performance in 
these asymptotic regions. Thus, under typical fading conditions, the 
transmitted carrier power might be increased at will without improving 
the ber significantly. The effects of decreasing the carrier power are 
discussed in Section 4.4. 

A family of curves like those shown in Fig. 9 was obtained (but are 
not given here) for each of nine uniformly spaced frequency offsets 
below midchannel to characterize the digital radio system sufficiently 
for the prediction of outage. A number of spot checks were also made 
using both rf and if fade simulators at symmetrical positive and 
negative offset frequencies, to establish that acceptable symmetry 
existed. 
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IV. CALCULATION OF OUTAGE 

This section describes four methods of calculating outage. The 
derivation of the critical curves of A and B, which provide the basis 
for making and understanding these calculations, is given in Section 
4.1. In Section 4.2 the detailed calculation of outage from the critical 
A-B curves is described. It is shown in Section 4.3 that for the present 
system this method may be greatly simplified by calculating only 
selectivity-caused outage (i.e., neglecting thermal noise). Section 4.4 
presents an approximate method of accounting for the effects of 
thermal noise. Section 4.5 provides a basis for estimating the selectiv- 
ity-caused outage from a single measurement. 

4. 1 Derivation of critical characteristics 

To calculate the outage for a fixed bit error rate, one must first 
obtain the critical curves of A and B at each simulated value of f , the 
notch position. Thus, from Fig. 9 which corresponds to /„ = -19.8 MHz 
(or <f> = -45°), we obtain six points on the critical curve of A and B for 
a ber of 10~ 3 , one point from each of the six curves which cross the 
critical ber. The value of B is obtained from the value of b since 

B = -20 log(l - b). (8) 

For the curve in Fig. 9 corresponding to B = 4.4 dB, we obtain the 
corresponding critical value of A for a ber of 10~ 3 from the value of 
carrier-to-noise ratio, which is 20.2 dB where this curve crosses the 
10 -3 ber line. Since the carrier-to-noise ratio is 62 dB when the channel 
is unfaded, the 20.2 dB value corresponds to a relative average power 
loss of 41.8 dB, 

L a = 62 - 20.2 = 41.8 dB. (9) 

Without loss of generality, we assume that the psk signal has a 
rectangular spectrum of width /&; consequently, the relative power 
transmitted by the model is obtained from eq. (2) as* 



i r fb 



I Ur-?,b co, 2^,r[^^)\. (10) 



* The calculated result is not critically dependent on the flatness of the signal 
spectrum or the spectral width chosen. We have used for f b a value of 25.3 MHz as 
representing the effective width of the signal. 
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Defining a correction term by 

C = -10 log jl + b 2 -2b cos 2tf oT fe-^ j, (11) 

we may express the signal loss as 

L s = -10\ogP ai , = A + C. (12) 

Thus, we obtain the critical value of A as 

A = L S - C. (13) 

For B = 4.4 dB (b = 0.4) and f = -19.8 MHz, we find C = 2.06 dB and 
the critical value of A is 41.8 - 2.1 = 39.7 dB. 

Carrying out these calculations for the six curves in Fig. 9, one can 
generate the critical curve of A and B for f = -19.8 MHz and a ber 
of 10~ 3 . The curve is shown in Fig. 10 along with the critical curves for 
several other values of the ber. A complete set of curves must be 
generated for all values of /„. 

The curves in Fig. 10 are typical of the critical curves obtained for 
| f 1 < 33 MHz. The intercept with the A-axis represents the flat fade 
margin for the given ber; this margin is independent of notch position. 
The intercept of a critical contour with the B-axis represents the 
shape, or relative fade depth, margin for the given notch position. For 
values of B to the right of this intercept, the critical value of ber 
cannot be obtained at any carrier-to-noise ratio for the given notch 
position. 

4.2 Outage calculation — detailed method 

The probability, P , of finding A and B outside all critical contours 
may be written with eqs. (5), (6), and (7) as 



■£ 



P*(<t>)Pc(4>) d<}>, (14) 

J — n 

where 



Pc(<l>)=\ P a(Y) P b(X) dY dX, (15) 

'0 JA C (X) 

and A C (X) is the functional relation of the critical values of A to B (or 
X), for B less than B c , the B-axis intercept, and for a given ber and <f> 
value.* Since measurements were made for a uniformly spaced set of 
notch positions with spacing A<> = 10°, we may approximate (14) by 



* The dependence of the function AAX) and the asymptote B, on ber and <f> is not 
explicitly denoted to keep notation simple. 
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Fig. 10— Critical curves of A and B for f„ = -19.8 MHz. 

All <*>, 



(16) 



To illustrate the calculation of outage probability with eqs. (15) and 
(16), we shall calculate the term in the summation of eq. (16) corre- 
sponding to a ber of 10" 3 and <f>, = -45° (or /„ = -19.8 MHz). From 
Fig. 11, which is taken from Fig. 10, we note that the double integral 
in eq. (15) may be broken into integrations over two regions. Thus 



PA*,)-* I J P a(Y) P b(X) dY dX 



(17) 



+ 



B, 



p A (Y)p B (X) dY dX, 



JAJXi 



where the two double integrals correspond to integrations over Regions 
1 and 2, respectively, in Fig. 11. Outage due to the occurrence of A and 
B in Region 1 may be described as outage due only to shape or 
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selectivity. In Region 2, outage is due to the combined effects of signal 
loss and selectivity. 

Using eqs. (5) and (6), the integral over Region 1 is obtained as 
e -B c /3.a rpj^ con tnbution due to thermal noise and shape (Region 2) is 
slightly more complicated. Dividing the interval to B c in Fig. 11 into 
N subintervals, as shown in Fig. 12, the probability of being in Region 
2 is the sum of the probabilities for each subinterval. Thus eq. (17) 
becomes 



Pcfo.) = e-^ /38 + £ [<r fl *- /38 - e- B * /38 ]P g ( Ak Ao{Bk) \ 



(18) 



where 



PAX) 



--f 



,-^/2 



dx. 



(19) 
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Fig. 1 1 — Classification of outage with respect to critical curve for ber = 10~ 3 . f = 
19.8 MHz. 
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Fig. 12— Outage calculation for an incremental interval. 

Evaluating the two components of eq. (18) from Fig. 11, we find 

P c (-45°) = 0.181 + 0.003 = 0.184. (20) 

This calculation was performed for 10 values of <f>, from —5 to —85 
degrees in 10-degree steps. Using these results in eq. (16) and multi- 
plying by two to account for positive values of <#>, which are assumed to 
contribute equally, we find the probability, P a for a ber of 10~ 3 as 

P = 0.0996. 

The expected outage for the data base period is, then, 

To = 5400 X 0.0996 = 538 seconds. (21) 

4.3 Outage calculation — selectivity only 

It is apparent from eq. (20) that most of the outage for the system 
under study is caused by selectivity, fades characterized by A and B 
values in Region 1. From eqs. (14) and (17), we may express P «, the 
probability of outage due to selectivity, as 
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Poh = pMpb(X) dX d$. (22) 

J -v JB C 

For the system studied for a ber of 10~ 3 , a finite B c is obtained only 
for | <£, | < 90 degrees. Hence, we may use eq. (7) to simplify (22)* 

From eq. (23) we see that the outage due only to selectivity depends 
on the relationship between B c , the asymptote of critical B values, and 
the notch angle or notch frequency. Figure 13 shows the relationship 
between B e and the notch frequency for four values of ber. It is 
apparent from eq. (22) that the outage probability is the probability of 
finding B and /„ values in the region above this curve. Such curves, 
therefore, provide a useful basis for evaluating the selectivity outage 
of a digital radio system. 

4.4 Outage calculation — approximate method 

For a radio system sensitive to both thermal noise and selectivity, 
the calculation of Section 4.3 is inadequate and that of Section 4.2 is 
unduly cumbersome. 

To illustrate a simpler, but generally applicable, method and at the 
same time to provide a useful incidental result, let us evaluate the 
effect of reducing the transmitted power by 10 dB. For the reduced 
power system, the carrier-to-noise ratio would be 52 dB for the unfaded 
channel, and the critical curves of A and B would be shifted by 10 dB. 
Figure 14 shows the critical curve of A and B for a 10 -3 ber and /„ = 

— 19.8 MHz with an overplot of the conditional mean of the distribution 
of A. The dotted curves on Fig. 14 represent 2-sigma intervals on 
either side of the mean. From the properties of the Gaussian distri- 
bution, one may determine that more than 95 percent of the values of 
A and B will lie between these two dotted curves. We designate as A m 
and B m the coordinates of the intersection of two curves: the critical A 

— B curve and the conditional mean curve. Then approximating the 
critical curve of A and B with a straight line segment tangent at (A m , 
B m ), with slope s, we may approximate the probability of outage by 
integrating the probabilities over the region to the right of the tangent 
line. Using eqs. (15) and (16), we obtain 



P„ = A</> £ ftta) p A (Y)p B (X)dY dX. (24) 

AU *i Jo JA m +s{B-B m ) 



* The factor of two is required in eq. (23) because the indicated summation corre- 
sponds to an integration only over negative notch frequencies (<J>, < 0). 
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Fig. 13— Asymptotic performance curves. Locus of values off,, and B that produce a 
fixed ber at high carrier-to-noise ratio. 




^ A + 2u ^.-»' 



10 15 20 

B = -20 Log (1-b) IN DECIBELS 



25 



30 



Fig. 14— Approximate outage calculation for 10 dB less transmitted signal. 
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Interchanging the order of integration and ignoring* the B dependence 
of A (B), this becomes 

P«-A* I p t (<t>i)e- B " /3 *e 0M *'*. (25) 

All <j>, 

Evaluating eq. (25) for a 10~ 3 ber and multiplying the result by 5400 
gives an outage estimate for the data base period of 602 seconds. 
Recalculating the total outage time at a 10~ 3 ber using the method of 
Section 4.2 [eqs. (16) and (18)] gives 636 seconds, which verifies the 
accuracy of the approximate method. The estimate of 636 seconds was 
calculated as an upper bound; the 602 seconds calculated using (25) 
tend to be a lower bound. We conclude that backing off transmitted 
power by 10 dB would increase the outage by about 12 percent (538 to 
602). 

4.5 A further simplification 

In this section, we show that the outage due to selectivity can be 
estimated approximately for a given ber from a determination of the 
in-band selectivity required (with the notch out-of-band) to produce 
that ber. Such a measurement may provide a useful approximation 
for any digital system using quadrature modulation components; 9 
however, we provide a justification based on the performance of the 
system at hand. In-band selectivity is defined as the difference between 
the maximum and minimum attenuation present in the (25.3-MHz) 
channel bandwidth. 

Since the in-band selectivity is a constant for any of the curves 
shown in Fig. 9, one can use Fig. 9 to plot the asymptotic ber against 
in-band selectivity for f Q = -19.8 MHz. Such a plot was generated for 
each notch position measured to produce the family of curves shown 
in Fig. 15. Note that, except for notch positions near the band center, 
the ber is uniquely related to the in-band selectivity. Neglecting the 
in-band notches, we find that a 10~ 3 ber corresponds to an in-band 
selectivity of 5.5 dB. 

If we use the original model of eq. (2) to determine the values of B 
that will produce an in-band selectivity of 5.5 dB for a number of 
different notch positions, we would generate Fig. 16. It is apparent that 
for this system there is a good correspondence between the curves of 
asymptotic performance (Fig. 13) and the curves of constant in-band 
selectivity (Fig. 16). 

To reinforce this conclusion, we provide Figs. 17, 18, and 19. Figure 
17 shows the locus of in-band selectivity in a 25.3-MHz band corre- 



* Including the effect of the slope of A (B) at B = B m gives the same symbolic result 
with s interpreted as the algebraic sum of the slope of the tangent and dAJdB evaluated 
at B = B m . 
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Fig. 15 — Measured asymptotic bit error rate vs peak-to-peak amplitude difference in 
a 25.3-MHz band. 




-10 -20 -30 -40 -50 

f Q , DISPLACEMENT OF NOTCH FROM MIDBAND IN MEGAHERTZ 



-60 



Fig. 16— Locus of B and /„ for modeled fades that have fixed peak-to-peak amplitude 
in a 25.3-MHz band. 



sponding to each of the curves of constant ber in Fig. 13. That is, for 
each ber and each value of notch position, f a , we have plotted the 
peak-to-peak amplitude difference in the band for the corresponding 
value of B c , the asymptotic critical value of B. Figure 18 shows a 
similar set of curves with the peak-to-peak delay distortion in a 25.3- 
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Fig- 17 — In -band selectivity (in 25.3-MHz bandwidth) corresponding to asymptotic 
critical values of notch depth (B c ) for several values of ber. 
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Fig. 18 — Peak-to-peak envelope delay distortion in 25.3-MHz bandwidth correspond- 
ing to asymptotic critical values of notch depth (B L ) for several values of ber. 



MHz band as the ordinate. Similarly, Fig. 19 has as the ordinate the 
"slope," or amplitude difference at a separation of 25.3 MHz. It is again 
clear from these three figures that the in-band selectivity is the 
relevant channel impairment giving rise to errors. We see from Fig. 18 
that, for out-of-band notches, high bers are obtained with very small 
values of peak-to-peak delay distortion, and from Fig. 19 that for in- 
band notches high bers are obtained for very small values (zero at 
mid- band) of slope. 
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Fig. 19— Amplitude difference at a 25.3-MHz frequency separation corresponding to 
asymptotic critical values of notch depth (£«•) for several values of BER. 

The model data base was analyzed to determine the time during 
which the in-band selectivity in a band of 25.3 MHz exceeded a given 
value. Figure 20 presents this distribution for in-band selectivity as 
calculated from the modeled fits. Figure 20 can be used directly in 
conjunction with Fig. 15 to calculate the outage times for the model 
data base.* For instance, from Fig. 15 we note that 5.5 dB of selectivity 
corresponds to a 10~ 3 ber. We use Fig. 20 to determine that 5.5 dB 
was exceeded for 520 seconds. 

V. COMPARISONS OF CALCULATED AND OBSERVED OUTAGES 

Using the methods of Sections 4.2 to 4.5, outage times were calcu- 
lated for bit error rates of 10" 3 to 10" 6 for both the model data base 
period and for a heavy fading month, by multiplying the outage 
probabilities by 5400 and 8100, respectively. 

5. 1 Model data base period 

Calculated and observedf outages for the model data base period 
are shown in Table I. In general, comparing the calculated results with 
observed results, we see that the outage is underestimated at high 
bers and overestimated at low bers. Any estimation procedure based 
on the current modeled state of the channel will tend to underestimate 



* In practice, one would use a single measurement of in-band selectivity. For instance, 
in Fig. 9 one would take the 5.7-dB value corresponding to the curve asymptotic at a 
10~ 3 BER. 

t Because of quantization, the outage times observed from the field experiment 
correspond to bit error rates of 1.26 X 10~ 3 , 1.57 x 1(T\ 0.981 X 10 r \ and 1.19 X NT 6 . 
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DECIBELS 

Fig. 20— Distribution of in-band selectivity (25.3-MHz bandwidth) for model data 
base. 

Table I — Outage in modeling data base period (seconds) 



10" 



io- 



10" 



10" 



Observed 

Detailed calculation (Section 4.2) 
Approximate calculation (Section 4.4) 
Asymptotic calculation (Section 4.3) 
Selectivity calculation (Section 4.5) 



636 


903 


1191 


1487 


538 


960 


1430 


1860 


527 


950 


1420 


1830 


527 


950 


1420 


1830 


510 


900 


1570 


2730 



outage at high bers because of hysteresis effects in the radio receiving 
equipment. That is, when the channel condition becomes sufficiently 
severe, the bit error rate becomes high enough (on the order of 10 -3 ) 
that the timing and/or phase of the radio system breaks lock. If the 
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channel impairment becomes less severe, the ber will not improve 
until the system resynchronizes. The hysteresis is important at the 
10 -3 ber, since a significant fraction of the events that cause 10 -3 ber 
will cause the system to break lock. 

One would expect to overestimate the outage at low bers because 
of the method of taking data. Recall that, in measuring the curves in 
Fig. 9, it was found that the ber depended on the phase to which the 
system had locked. The recorded performance represented the worst- 
phase condition. At a 10" 6 ber, the best phase produces a ber that is 
about Vz that produced by the worst phase; the difference in ber from 
worst to best phase at a 10 -3 ber is negligible. Hence, one would 
expect outage to be overestimated significantly at low bit error rates. 

In comparing the outage calculated from in-band selectivity (Section 
4.5) to the outage observed, we find that the overestimation of outage 
at low bers is more severe than with the other methods. This is due 
to the greater sensitivity of the differential selectivity method to the 
bias induced by choosing the worst-case phase. For instance, compar- 
ing calculations at a 10" 6 ber, we find that Fig. 20 is steeper for 
amplitude differences near 2 dB than is Fig. 3 near B values of 3.5 dB. 
(Figure 9 verifies the appropriateness of this comparison). More gen- 
erally, one expects the method based on in-band selectivity to overes- 
timate the outage because the method is based on notches out of band. 
From Fig. 15, it is apparent that, for a given AA, some scans will not 
have the ber specified. 

We conclude that, although calculation of outage from sensitivity to 
in-band selectivity provides quick estimates, they are less accurate. 
The calculation requires knowledge of the distribution of in-band 
selectivity over a specified bandwidth. These statistics are neither 
simple nor generally available. It has been shown, 10 for instance, that 
slope statistics have a nontrivial dependence upon the measurement 
bandwidth. 

It is clear that the calculations based on selectivity (Sections 4.3 and 
4.5) agree for the system studied here because that system has very 
little outage due to thermal noise limitations, and because it is sensitive 
primarily to in-band amplitude excursions. The extent to which these 
statements are true for other systems is currently unknown. 

5.2 Outage on a monthly basis 

The results in Table I may be put on the basis of a heavy fading 
month by increasing them by a factor of 1.5, as discussed in Section 
2.2. The resulting outages (including the scaled observed outage) are 
compared with the outage observed in a one-month period 8 in Table 
II. We see that the outage times observed in the total one-month 
period agree well with the values obtained by scaling the observed 
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Table II — Outage in 


a heavy fading month (seconds) 






BER = 10" 3 


io- 


10 -5 


lO -8 


Observed (Ref. 8) 
Scaled observation from Table I 
Calculation (Section 4.2-4.4) 
Selectivity calculation (Section 4.5) 


1000 
955 
800 
770 


1320 
1350 
1430 
1350 


2100 
1790 
2140 
2350 


2900 
2230 
2760 
4100 



outage for the data base period used in modeling, except for the slight 
divergence appearing at low bers. This divergence should not be 
unexpected for this equipment. As may be seen in Fig. 15, a 10~ 6 ber 
may be caused by differential amplitude selectivity in band of 2 dB. 
Such modest amounts of selectivity may be expected to occur some- 
times in the presence of very moderate selective fading. The modeling 
data base was constructed by selecting only periods of significant 
selective fading. This reinforces the comments made in conjunction 
with Fig. 3, namely, that the model distribution of B represents a lower 
bound for small values of B which can contribute to outage at the 10" 6 
ber level. 

VI. CONCLUSIONS 

We have demonstrated the validity of a technique for estimating the 
unprotected outage of a digital radio system due to selective fading on 
a particular hop in the 6-GHz common carrier band. The technique 
required field measurements to statistically characterize the param- 
eters of a model of propagation on the hop. It also requires performance 
data obtained in the laboratory for the radio system by stressing it 
with a two-path fade simulator with a differential path delay of 6.3 ns, 
corresponding to the fixed delay channel model. Since the radio path 
on which these measurements were made has a length close to the 
average for the Bell System long haul radio network and has an 
average incidence of fading activity, the channel model is representa- 
tive of a typical path. At the very least, the technique provides a basis 
for determining the relative merits of various digital radio systems 
operating without benefit of space diversity. 

For the system under test, outage was calculated by four different 
methods. Because this system was selectivity-limited rather than 
noise-limited, all four methods predicted approximately the same 
outage as that summarized in Table I; however, the method based on 
in-band selectivity is more severely biased at low bers. The method 
based on asymptotic performance and that based on in-band selectivity 
can only be used to estimate outage due to selectivity. If the transmit- 
ted power of the system under test were reduced by 10 dB, both of the 
other two methods, the detailed and the approximate method, would 
predict an increase in outage time of about 12 percent. 
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